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ABSTRACT 

Spectral lag, the time difference between the arrival of high-energy and low-energy photons, 
is a common feature in Gamma-ray Bursts (GRBs). Norris et al. 2000 reported a correlation 
between the spectral lag and the isotropic peak luminosity of GRBs based on a limited sample. 
More recently, a number of authors have provided further support for this correlation using 
arbitrary energy bands of various instruments. In this paper we report on a systematic extraction 
of spectral lags based on the largest Swift sample to date of 31 GRBs with measured redshifts. 
We extracted the spectral lags for all combinations of the standard Swift hard x-ray energy 
bands: 15-25 keV, 25-50 keV, 50-100 keV and 100-200 keV and plotted the time dilation corrected 
lag as a function of isotropic peak luminosity. The mean value of the correlation coefficient for 
various channel combinations is -0.68 with a chance probability of ~ 0.7 x 10~ 3 . In addition, the 
mean value of the power-law index is 1.4 ± 0.3. Hence, our study lends support for the existence 
of a lag-luminosity correlation, albeit with large scatter. 

Subject headings: GRB; gamma rays bursts; redshift, spectral lag 



1. Introduction 

After decades of research, a satisfactory expla- 
nation of the temporal behavior of Gamma-ray 
Burst (GRB) light-curves is still lacking. Despite 
the diversity of GRBs, some general characteristics 
and correlations have been identified: Spectral lag 
is one such characteristic. The spectral lag is the 
difference in time of arrival of high-energy pulses 
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verses low-energy pulses. In our analysis a positive 
spectral lag corresponds to an earlier arrival time 
for the higher energy photons. The o bserved spec- 
tral lag is a common feature in GRBs ( Cheng et al.l 
19951: iNorris et al]|l996l: iBand Ill997h . The study 
of spectral lag between energy bands, which com- 
bines temporal and spectral infor mation, poten- 
tially can constr a in GRB models (|Lu et al.l [2006: 
Shen et al.l [20051 iQin et al.l r2003 Ischaeferl 12004 ; 



Ioka fc Nakamurall200lt ISalmonson Il2000h 

Based on six G RBs with known redshifts, 
Norris et al. ] (I2000T I found an anti-correlation be- 
tween the spectral lag and the isotropic peak lu- 
minosity. Furthe r eviden c e for t his correlation 
pr ovided b y Norris ( 20021). Gehrels et al. 
6l).ISchaeferl (120071) . IStamatikos et al.l (|2008aF ) 



(2006 



and Hakkila et al 



(|200c€ 



Others have used this 



relation as a redshift indicator (Murakami et al 



2003; Band et al. 


20041) and as a cosmological 


tool dBloom et al. 


20031: ISchaeferl2007: Liang et al. 



2008; Mos quera Cuesta et al.l 



Hak kila et al. I (120081) have used a pulse-profile 



fitting technique (a four-parameter pulse function 
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introduced by iNorris et al. ( 20051) 1 to show that 
the correlation is between lags of the pulses and 
the luminosity of the pulses seen in GRBs. How- 
ever, the method is limited because it applies only 
to very bright bursts where pulses are clearly iden- 
tifiable and described by the assumed pulse profile. 

Many authors have tried to explain the physical 
cause of the lag-luminosity rel ation and a number 
of models have been proposed. ISalmonson I d2000h 



argues that the anti-correlation is due to the varia- 
tions in the line-of - sight velocity of various GRBs. 
Ioka fc Nakamural (|200ll ) suggest that the relation 
is a result of variations of the of f - axis a ngle when 
viewing a narrow jet. ISchaeferl fj2004h invokes a 



rapid radiation cooling effect to explain the corre- 
lation. This effect tends to produce short spectral 
lags for highly luminous GRBs. 

Regardless of its physical origin, the spectral 
lag is an important measurement for GRB science 
because of its us efulness in differentiatin g long 
and short GRBs (|Kouveliotou et al.lll993l ): Long 
bursts give large lags and shor t bursts give rela- 
tively small lags ( Norris 19951: Norris fc Bonnell 
2006). Even though a few exceptions to this clas- 
sification scheme have been found, such as GRB 



060614 (jGehrels et al.ll2006l ). the GRB community 
still continues to use the spectral lag as one of the 
classification criteria. Note that more elaborate 
classification schemes based on multiple observa- 
tional parameters, such as t he host galaxy prop- 
erty, has also been proposed (jDonaghv et al. 2006; 



Zhang et alj |2009) 



Moreo ver, based on the an alysis of GRB 
080319B. IStamatikos etafl (|2009l) show that there 
is a possible correlation between the prompt opti- 
cal emission and the evolution of spectral lag with 
time. 

Most of the previous work on spectral lags has 
been based on observations with the Burst and 
Transient Source Experiment (B ATSE) on the 
Compton Gamma Ray Observatory dTsutsui et al 
«T It 



2008t lHakkila et al J [20081 120071: [Chen et al.l | 2005 
Band et al.l |2004 ISalmonson fc Galamal 12002 : 
Norris I l2002t iBand I Il997l). The launch of the 
Swift satellite (jGehrels et al.l 120041) ushered in a 



new era of GRB research. In this paper we present 
a detailed study of spectral lags using a subset of 
Swift Burst Alert Telescope (BAT) data. 

The structure of the paper is the following: In 



section [2] we discuss our methodology with a case 
study featuring GRB 060206. In section [3] we 
present our results for a sample of 31 Swift BAT 
long bursts and investigate the lag-luminosity re- 
lation for various channel combinations. Finally, 
in section [4] we discuss some implications of our 
results. Throughout this paper, the quoted uncer- 
tainties are at the 68% confidence level. 

2. Methodology 

2.1. Light Curve Extraction 

Swift BAT is a highly sensitive instrument, 
which utilizes a coded apertu re mask to localize 
bursts (jBarthelmv et al.1l2005 ) . The basic imaging 
scheme is that a gamma-ray source illuminates the 
coded aperture mask, and casts a shadow onto a 
position sensitive detector. Each position in the 
sky will produce a unique shadow pattern in the 
detector plane. Hence by comparing the observed 
shadow with precalculated shadow patterns for all 
possible points in the sky it is possible to find the 
actual position of the source that created the given 
shadow pattern. However, in practice each detec- 
tor can be illuminated by many points on the sky 
whereas each point on the sky can illuminate many 
detectors. To disentangle each point in the sky, 
special software designed by the Swift BAT team 
is used. 

In order to generate light curves, a process 
called mask weighting is utilized. The mask 
weighting assigns a ray-traced shadow value for 
each individual event, which then enables the 
user to calculate light curves or spectra. We 
used the batmaskwtevt and batbinevt tasks in 
FTOOLS to generate mask weighted, background- 
subtracted light curves for our analysis. Resulting 
light curves and their uncertainties are calculated 
by propagation of errors from raw counts (subject 
to Poissonian noise). 

2.2. The Cross Correlation Function and 
Spectral Lag 

There are at least three well known ways 
of extra cting spectral lag s : (1) pulse peak-fi t 



method (Norris et al. 2005t Hakkila et al 



2008), 



(2) Fourier analysis method (|Li et al.l 120041 ). and 

(3) cross -correlation functi o n (CCF) a nalysis 
method (|Cheng et al.l Il995t IBand I Il997h . The 
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pulse peak-fit method gives a simple straight for- 
ward way for extracting lags. It does however 
assume a certain pulse function for the pulses in 
the light curve and may also be limited to very 
bright bursts. It is not immediately clear how 
this method would fare in cases where the light 
curves are sufficiently complex i.e., not dominated 
by a prominent pulse. For transient events such 
as GRBs using the Fourier analysis technique also 



has its difficulties ( Li et all 120041 ) . Since GRB 



light curves do not exhibit obvious periodicities, 
Fourier transforms typically yield a large number 
of coefficients to describe their temporal structure. 
These coefficients, in turn, produce a spectral lag 
value for each corresponding frequency component 
i.e., a spectrum of lags is generated. The gener- 
ated spectra exhibit a variety of shapes depend 



ing o n the complexity of the light curve (|Li et al 



2004) thus making the extraction of an intrinsic 
lag questionable. Hence, in this work, we develop 
a method to calculate the time-averaged spectral 
lag and its uncertainty via a modification of the 
CCF method. 

The use of the Pearson cross-correlation func- 
tion is a standard method of estimating the degree 
to which two series are correlated. For two count- 
ing series Xi and xji where i = 0, 1, 2, ...(N — l), the 
CCF with a delay d is defined as 



CCF std (d,x,y) 



x)(yz+d - v) 



(i) 

Here x and y are average counts of the two series 
x and y respectively. The denominator in the ex- 
pression above serves to normalize the correlation 
coefficient such that —1 < CCF std (d, x,y) < 1, 
the bounds indicating maximum correlation and 
zero indicating no correlation. A high negative 
correlation indicates a high correlation but of the 
inverse of one of the series. Note that the time 
delay (r) is given by r = d x time bin size. 



However. I Band I ([19971 ) proposed that for tran- 
sient events such as GRBs, non-mean subtracted 
definition given below is more suitable for the 
time-averaged lag. 



CCF Bajld (d,x,y) 



Emm(N.N-d) 
i=max(l,l-d) Xi Vi+d 



(2) 



VE* x'i E l y't 

We have tested both definitions of the CCF us- 
ing synthetic light curves with artificially intro- 



duced spectral lags. Our tests showed that the 
CCF Band consistently recovered the introduced lag 
while CCF std sometimes f ailed (possibl e reasons 
for this failure are noted in iBand I (|l997l )). Hence 
in our analysis we used the CCF Band definition and 
from this point onwards in the paper we refer to 
it simply as the CCF. 

For a given pair of real light curves, we deter- 
mine the CCF using Equation[2J At this stage the 
resulting CCF values do not have any uncertain- 
ties associated with them. In order to determine 
these uncertainties, we use a Monte Carlo simu- 
lation. Here we make 1,000 Monte Carlo realiza- 
tions of the real light curve-pair based on their 
error bars as shown below. 



LC 



simulated 
bin 



= LC 



real 
bin 



CxLcg 



cal c 



(3) 



Where £ is a random number generated from a 
gaussian distribution with the mean equal to zero 
and the standard deviation equal to one. For each 
simulated light curve-pair we calculate the CCF 
value for a series of time delays. This results in a 
1,000 CCF values per time-delay bin. The stan- 
dard deviation of these values per time-delay bin 
is then assigned as the uncertainty in the original 
CCF values obtained from the real light curves. 

2.3. Extracting Spectral Lags 

We realize there may be number of ways to de- 
fine the spectral lag, but in this work, we define 
it as the time delay corresponding to the global 
maximum of the cross-correlation function. To lo- 
cate this global maximum, we fit a Gaussian curve 
to the CCF. The uncertainties in the CCF are ob- 
tained using a Monte Carlo procedure discussed 
in section 12.21 In essence, our fitting procedure 
locates the centroid of the cross correlation func- 
tion and is thus relatively insensitive to spurious 
spikes in the CCF. We tested and verified the ro- 
bustness of this procedure by performing a num- 
ber of simulations in which artificial lags were first 
introduced into the light curves and then success- 
fully recovered. In addition, our tests with these 
artificial light curves show that the CCF can be- 
come asymmetric (around its global maximum) if 
the shape of one of the light curves is significantly 
different from the other. This energy dependent 
feature potentially requires a more complex fitting 
function than a Gaussian or a quadratic to fit the 
CCF over the entire range. Instead of resorting 
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to a more complex fitting function we were able to 
recover the (known) lags by fitting the CCFs (with 
a Gaussian) over limited but asymmetric ranges. 

2.3.1. Time Bin Selection 

For Swift GRBs the minimum time binning 
is 0.1 ms but one can arbitrarily increase this all 
the way up to the duration of the burst. It is 
important to understand the effect of time binning 
on the extracted spectral lags. Presumably, by 
changing the time binning of the light curve one is 
affecting the signal-to-noise ratio. By employing 
increasingly coarser binning one is averaging over 
the high frequency components of the light curve. 
Clearly, one has to be careful not to use overly 
large time bin sizes otherwise one risks losing the 
sought-after information from the light curve. 

In order to understand the effect of time bin- 
ning more fully, we did a number of simulations 
utilizing peak normalized synthetic light curves 
(composed of FREE0-like pulse shapes with Gaus- 
sian distributed noise) in which artificial lags were 
introduced. We incrementally increased the noise 
level and studied its effect on the maximum cor- 
relation value in the CCF vs time delay (CCF- 
Max) plot. In Fig. Q]we display the synthetic light 
curves with several noise levels (0%, 20%, 40% and 
60% respectively) as well as the calculated CCF 
with typical Gaussian fits. As expected, the CCF- 
Max value (see the right panel of Fig. [TJ decreases 
gradually as the noise level increases. We also note 
that the scatter in the CCF increases considerably 
with the noise level. The global maximum in the 
CCF is clearly visible at the 40% noise level and 
a good fit is obtained with a Gaussian. However, 
this is not the case for the 60% noise-level curve, 
in which the scatter is quite significant, and the 
CCF global maximum is barely visible leading to 
a poor fit. 

In Fig. [5] we show the behavior of the CCFMax 
value and the extracted spectral lag as a function 
of the noise level. We first note that the CCFMax 
value smoothly tracks the signal-to-noise level in 
the light curves (see the upper panel of Fig. [2]). 
Secondly, we note that the extracted lag value 
agrees well with the artificially introduced lag of 
10 seconds up to a noise level of about 40%. We 
further note that the scatter in the extracted lag 



value increases as the noise contribution increases 
beyond 40%. Although it is not immediately ob- 
vious from this figure, CCF values above a noise 
level of 40% show large scatter (see bottom right 
panel of Fig. [T} thus making the extracted lag 
value uncertain. This is directly reflected by the 
increasing error bars in the extracted value. These 
simulations were repeated for number of time lags 
and in all cases similar results were obtained, in 
particular, the behavior of the CCFMax as a func- 
tion of the noise level was confirmed. Based on the 
results of these simulations, we chose a CCFMax 
~ 0.5, corresponding to a noise level of about 40%, 
as our guide for picking the appropriate time bin- 
ning. 

Procedurally, we start with a time bin size of 
1024 ms and decrease the time binning by powers 
of two until the CCFMax becomes ~ 0.5 and use 
that time bin size as the preferred time binning 
for the lag extraction. By using this procedure 
we are able to arrive at a reasonable bin size that 
preserves the fine structure in the light curve and 
at the same time keeps the contribution of the 
noise component at a manageable level. 

2.3.2. Uncertainty in Spectral Lags 

We have studied three methods to determine 
the uncertainty in the extracted spectral lags. The 
first method is to use the uncertainty that is ob- 
tained by fitting the CCF with a Gaussian curve. 
The secon d method is an adaptation of equation 
(4) used in lGaskell fc Peterson I (jl987l) 



CTlag 



0.75 Whwhm 



1 + hy/n^2 



(4) 



iFast Rise Exponential Decay 



Here Whwhm is the half-width half-maximum of 
the fitted Gaussian, h the maximum height of the 
Gaussian and n is the number of bins in the CCF 
vs time delay plot. This method utilizes more in- 
formation about the fit and the CCF such as the 
width, height, and number of bins to estimate the 
uncertainty. The third method utilizes a Monte 
Carlo simulation. We found that the first method 
gives systematically smaller uncertainty in the lag 
by a factor of two or more relative to the other 
two methods. The second and third methods give 
comparable values. We adopted the most conser- 
vative of the three methods (i.e. the one based 
on the Monte Carlo simulation) to determine the 
uncertainties in the lag. 
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2.3.3. Lag Extraction: Case Study 

To illustrate the lag extraction procedure more 
clearly, we present a case study using GRB 060206. 
The light curve segment is selected by scanning 
both forward and backward directions from the 
peak location until the count rate drops to less 
than 5% of the peak count rate (using 15 — 200 
keV light curve). This selection method is cho- 
sen to include the most intense segment of the 
burst and to capture any additional overlapping 
pulses near the main structure. Presumably, these 
pulses also contribute to the overall spectral lag. 
In the case of GRB 060206 this corresponds to a 
light curve segment starting 1.29 seconds prior to 
the trigger and 8.18 seconds after the trigger (see 
Fig. [3]). Next we calculate the CCF and plot it 
as a function of time delay as shown Fig. |4] Error 
bars on the CCF points were obtained via a Monte 
Carlo simulation of 1,000 realizations of the orig- 
inal light curves (see section |2T2"|) . As noted ear- 
lier, we start with time bin size of 1024 ms and 
decrease the time binning by powers of two until 
the CCFMax becomes ~ 0.5 for a given channel 
combination, in this case BAT standard channel 2 
(25-50 keV) and 3 (50-100 keV). For GRB060206 
channel 2 and 3, the time bin size correspond to 
8 ms. The global maximum of the CCF vs time 
delay plot corresponds to the spectral lag and its 
value is obtained by fitting a Gaussian curve. We 
choose a range of the time delay (in this case from 
-1.5 seconds to 1.5 seconds) manually to identify 
the global maximum. In order to obtain the un- 
certainty in the spectral lag, we employ another 
Monte Carlo simulation, in which we create 1,000 
additional realizations of the input light curves as 
described in section [2"T2l and repeat the previously 
described process for the simulated light curves. 
A histogram of the resulting (1,000) spectral lag 
values is shown in Fig. for GRB 060206. The 
standard deviation of these values is the uncer- 
tainty in the spectral lag. 

2.4. Isotropic Peak Luminosity 

To compare observations with different instru- 
ments we need to calculate flux over some fixed 
energy band. In order to do this we need to 
know the best-fit spectral function to the ob- 
served spectrum and its spectral parameters. Of- 
ten, GRB spectra can be well fitted with the Band 



function (jBand et al.lll9 93"). an empirical spectral 
model defined as follows: 



TV = 



I U2+aH00kcVJ B ' ° lLX 



E n 



100 kcv 1 L(2+a)100keV 



rwise. 



There arc four model parameters in the Band func- 
tion; the amplitude (A), the low-energy spectral 
index (a), the high-energy spectral index (J3) and 
the peak of vF v spectrum (E p ). 

If the GRB spectrum is well described by the 
Band function, then the values of a, /?, E p and 
the observed peak flux, / bs, in a given energy 
band (£mi n and E max ) are often reported. We 
can calculate the normalization A with 



obs 



j*-;N>(E)dE 



(5) 



where / Q b s is given in photons cm 2 s 1 and N' = 
N/A. 

The observed peak flux for the source- frame en- 
ergy range E x = 1.0 keV to E 2 = 10, 000 keV is 



r E 2 /(l+z) 

/X w =/ N(E)EdE. 

The isotropic peak luminosity is 

L iso = 47rd|/X w 
where dh is the luminosity distance given by, 



d L = 



(l + z)c 



dz' 



Ho Jo y/n M (i + z'T + ni 



(6) 



(7) 



(8) 



For the current universe we have assumed, Qjvf = 
0.27, Ol = 0.73 and the Hubble co nstant Hp is 



70(kms~ 1 )/Mpc = 2.268x 10" 18 s" 1 ijKomatsu et al 



2009) 



To determine the uncertainty in Li so , we employ 
a Monte Carlo simulation. We simulate spectral 
parameters a, j3, E p and flux assuming their re- 
ported value as sample mean and reported uncer- 
tainty as sample standard deviation, then calcu- 
late iiso, for 1,000 variations in these parameters. 
If a parameter has uneven uncertainty values then 
each side around the parameter is simulated with 
different uncertainty values as standard deviation. 
Then we take the 16th and the 84th ranked values 
(la uncertainty) as the lower limit and the upper 
limit of Lj so respectively. 
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3. Results 

We selected a sample of long GRBs (T90 > 2 
sec, excluding short bursts with extended emis- 
sion), detected by Swift BAT from 2004 Decem- 
ber 19 to 2009 July 19, for which spectroscopi- 
cally confirmed redshifts were available. Out of 
this initial sample (102), a subset of 41 GRBs 
were selected with peak rate > 0.3 counts/sec/det 
(15 — 200 keV, 256 ms time resolution). Finally, 
we selected 31 GRBs for which a clear global max- 
imum can be seen in the CCF vs time delay plots 
with maximum correlation of at least 0.5 (with 
256 ms time binning) for all channel combinations. 
The spectral parameters of the final sample are 
given in Table [TJ We note that our final sample 
contains bursts with redshifts ranging from 0.346 
(GRB 061021) to 5.464 (GRB 060927) and the av- 
erage redshift of the sample is ~2.0. 

Out of our sample, 18 bursts have all Band 
spectral parameters measured and comprise our 
"Gold" sample. The remaining 13 bursts are fur- 
ther divided into "Silver" and "Bronze" samples. 
In the "Silver" sample, 10 bursts have E p deter- 
mined by fitting a cutoff power- law0 (CPL) to 
spectra and for GRB 060418, E p is reported with- 
out uncertainty, so we assumed a value of 10%. 
These 10 bursts do not have the high-energy spec- 
tral index, /3, measured, so we used the mean value 
of the BATSE dist r ibution, which is -2.36 ±0.31 
(jKaneko et al.l 120061: [Sakamoto et al"1l2009l ). The 
"Bronze" sample (consisting of 3 bursts) does not 
have a measured E p . We have estimated it us- 
ing the power-law index ( r) of a simple power-law 
(PL) fit as described in ISakamoto et al. ( 2009h . 
For these 3 bursts, the low-energy spectral in- 
dex, a and the high-energy spectral index, /3, were 
not known, so we used the mean value of the 
BATSE a and j3 distribution, which is -0.87±0.33 
and —2.36 ±0.31 resp ectively ( Kaneko et al. 20061 : 
Sakamoto et al.|[2009h . All estimated spectral pa- 
rameters are given in square brackets in Table [TJ 

Using the spectral parameters and redshift in- 
formation in Table [TJ we have calculated the peak 
isotropic luminosities for all the bursts in our sam- 
ple: these results are shown in Table [U GRB 
080430 has the lowest luminosity in the sample 
(~ 1.03 x lO^ergs- 1 ), and GRB 080607 has the 



highest luminosity (~ 7.19 x 10 53 ergs _1 ). The 
sample spans roughly three orders of magnitude 
in luminosity. 

We extracted the spectral lags for all combina- 
tions of the canonical BAT energy bands: channel 
1 (15-25 keV), 2 (25-50 keV), 3 (50-100 keV) and 
4 (100-200 keV). We took the upper-boundary of 
channel 4 to be 200 keV because we found that 
after the mask weighting the contribution to the 
light curve from energies greater than ^200 keV 
is negligible. The nomenclature is straightforward, 
i.e. the spectral lag between energy channels 4 and 
1 is represented by Lag 41. As such there are six 
channel combinations and the results for all six are 
shown in Table [3J The segment of the light curve 
used for the lag extraction (T + X§ and T + Xe, T 
is the trigger time) , the time binning of the light 
curve, and the Gaussian curve fitting range of the 
CCF vs time delay plot (with start time, and end 
time denoted as LS and LE respectively) are also 
given in Table [3] 

poooh . 



We noticed, as did Wu fc Fenimore 



that the lag extraction is sensitive to a number 
of parameters. Hence, in Table El we specify the 
band pass that we used to extract the lag, seg- 
ment of the light curve used, temporal bin reso- 
lution, and the fitting range used in the CCF vs 
time delay plot. These additional parameters are 
reported in order to facilitate reproduction of the 
results and direct comparison with other extrac- 
tion techniques. 

Figures [6] through [TTJ show log-log plots of 
isotropic peak luminosity vs redshift corrected 
spectral lag for various energy channel combi- 
nations. Red circles represent bursts from the 
"Gold" sample, blue diamonds shows bursts from 
the "Silver" sample and green triangles are bursts 
from the "Bronze" sample. The best-fit power- 
law curve is also shown in these plots with a dash 
line. Since there is a large scatter in these plots, 
to compensate, the uncertainties of the fit param- 
eters are multiplied by a factor of y/x 2 /ndi (see 
Table [4]) H. The dotted lines indicate the estimated 
1 a confidence level, which is obtained from the cu- 
mulative fraction of the residual distribution taken 
from 16% to 84%. 

It is interesting to note that GRB 080603B ex- 
hibits five negative lags out of six possible combi- 



2 dN/dE ~ E a exp (-(2 + a)E/E p ) 



3 ndf - number of degrees of freedom 
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nations. While these negative lags are not shown 
in the plots, it is worth noting that negat ive lags 
are not necessarily unphysical (|Rvdell2005l ) . More- 
over, in the few cases where the uncertainty is 
large, i.e. the extracted lags are consistent with 
zero, these points are not plotted either but are 
listed in Table [3] We recognize that the omission 
of these negative and zero lags is a potential source 
of bias. 

As seen in Fig. [5] through [TTJ our results sup- 
port the existence of the l ag- luminosity c orrela - 
tion originally proposed by iNorris et al. ] llxl). 
Table |4] lists the correlation coefficients for all 
six channel combinations. The lag for channel 
combination 31 has the lowest correlation with 
L; so , where the correlation coefficient is -0.60 (with 
chance probability of ~ 1.5 x 10 -3 ) and the lag 
for channel 43 has the highest correlation with 
coefficient of -0.77 (with chance probability of 
~ 3.0 x 10~ 4 ). However, we note that there is con- 
siderable scatter in the plots. The results of our 
best fit curves for each energy band combinations 
are also given in Table |4] The mean value of the 
power-law indices that we get for various channel 
combinations is 1.4 ± 0.3. O ur value is consistent 
with the 1.14 power-law index PNorris et al. ( 2000l) 



reported using lags between BATSE energy bands 
100 - 300 keV and 25 - 50 keV. Ou r result s are 
also cons i stent with IStamatikos et al.l (|2008a) and 
ISchaeferl (|2007l ) who reported values of 1.16 ± 0.21 
and 1.01 ± 0.10 (assuming an uncertainty of 10%), 
respectively. 

4. Discussion 



Band I (jl997h showed that gamma-ray burst 



spectra typically undergo hard-to-soft peak evo- 
lution, i.e. the burst peak moves to later times 
for lower energy bands. In our sample we have 
six lag extractions for each burst. The perfect 
hard-to-soft peak evolution scenario is indicated 
by positive lag values for all channel combinations 
plus lag41 > lag42 > lag43 and lag31 > lag32. 
However, all bursts i n our s a mple do not show 
this perfect behavior. iBand I (|l997l ) used a scor- 
ing method to quantify the degree of hard-to-soft 
peak evolution. We used a more elaborate scor- 
ing method to assign a score to each GRB as fol- 
lows: First, we increase the burst score by one if 
one of the six lag values is positive or decrease 



it by one if it is negative. Thus, a GRB can get 
a score ranging from -6 to +6 at this first step. 
Then, we compare the lag values of channel 4 as 
the base (lag43, lag42, and lag41). The score is 
increased by one if the burst meets one of the fol- 
lowing conditions: lag41 > lag42, lag41 > lag43 
or lag42 > lag43. We continue this procedure for 
channel 3 as the base also (lag31 > lag32). We de- 
crease the score by one if it is otherwise. According 
to this scoring scheme a score of +10 corresponds 
to the perfect case that we mentioned earlier. A 
positive score indicates overall hard-to-soft peak 
evolution in the burst to some degree. A negative 
value indicates soft-to-hard peak evolution. Out 
of 31 bursts in our sample 19 bursts show perfect 
hard-to-soft peak evolution with a score of +10. 
About 97% of bursts in our sample have a score 
of greater than zero, which is consis tent with the 
90% value reported bv lBandl (|l997h . 

If one wants to use the lag-luminosity relation 
as a probe into the physics of GRBs (in the source 
rest frame) , then a few corrections to the spectral 
lag are required; 1) correct for the time dilation 
effect (z-correction), 2) take into account the fact 
that for GRBs with various redshifts, observed en- 
ergy bands correspond to different e nergy bands at 
the G RB rest frame (k-correction) . iGehrels et al 



(|20061 ) approximately corrected observed spectral 
lag for the above mentioned effects. We also ex- 
amined these corrections. The z-correction is done 
by multiplying the lag value by (1 + z) . The 
k-correction is approximately done by multiply 
ing the lag value by ( 1 + z) 33 (IGehrels et al 



20061 IZhang et af1l2009h . In table [5] we list the 



correlation coefficients with no correction, only z- 
correction, only k-correction and both corrections 
applied. For example correlation coefficient of 
Lag31 and L- lso is -0.38 without any corrections. 
After the k-correction the correlation coefficient 
is -0.29. Therefore, we do not gain a significant 
improvement in the correlation by applying the 
k-correction. However, the correlation improves 
significantly after the z-cor rection (-0.60). The ap- 



proximate k-correction of ( Gehrels et al. 20061 ) 



is 



based on the assumption that the spectral lag is 
proportional to the pulse widt h and pulse width - 
is proportional to the energy (| Zhang et al. 2009: 
Fenimore et al.lll995 ). These approximations de- 



pend on clearly identifying a pulse in the light 
curve and may be of limited validity for multi- 
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pulse structures. A better method would be to de- 
fine two energy bands in the GRB rest frame and 
project those two bands into the observer frame 
and extract lags between them (Ukwatta, et al., 
20f0 in preparation). 

5. Conclusion 

In this work we have used the CCF technique 
to extract spectral lags for a sample of Swift 
BAT GRBs with known redshifts. By using Monte 
Carlo simulations, we have extended this tech- 
nique to reliably determine the uncertainties in 
the extracted spectral lags. Normally these un- 
certainties would be very difficult to calculate an- 
alytically. 

This study provides further support for the ex- 
istence of the lag- luminosit y corre lation, originally 
proposed by iNorris et alH (|2000h . We note how- 
ever, that there is a significant scatter in the cor- 
relation. 

The authors are indebted to the late Dr. David 
L. Band for fruitful and insightful discussions on 
the CCF methodology. In addition, we take this 
opportunity to acknowledge useful input from 
David A. Kahn regarding the luminosity cal- 
culations. We also thank the anonymous ref- 
eree for comments and suggestions that signif- 
icantly improved the paper. The NASA grant 
NNX08AR44A provided partial support for this 
work and is gratefully acknowledged. 
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Table 1 

grb redshift and spectral information 



GRB 


z 


Peak Flux a 


E P b 


a c 


/3 d 


Reference 



GRB050401 

GRB050603 

GRB050922C 

GRB051111 

GRB060206 

GRB060210 

GRB060418 

GRB060904B 

GRB060908 

GRB060927 

GRB061007 

GRB061021 

GRB061121 

GRB070306 

GRB071010B 

GRB071020 

GRB080319B 

GRB080319C 

GRB080411 

GRB080413A 

GRB080413B 

GRB080430 

GRB080603B 

GRB080605 

GRB080607 

GRB080721 

GRB080916A 

GRB081222 

GRB090424 

GRB090618 

GRB090715B 



2.899 1 
2.821 2 
2.199 3 
1.550 4 
4.056 5 
3.913 6 
1.490 7 
0.703 s 
1.884 9 
5.464 10 
1.262 11 
0.346 12 
1.315 13 
1.496 14 
0.947 15 
2.145 16 
0.937 17 
1.949 18 
1.030 19 
2.433 20 
1.101 21 
0.767 22 
2.689 23 
1.640 24 
3.036 25 
2.591 26 
0.689 27 
2.770 28 
0.544 29 
0.540 30 
3.000 31 



10.70 ± 0.58 
21.50 ± 0.67 
7.26 ± 0.20 
2.66 ± 0.13 

2.79 ± 0.11 
2.72 ± 0.18 
6.52 ± 0.22 
2.44 ± 0.13 
3.03 ± 0.16 
2.70 ± 0.11 

14.60 ± 0.23 
6.11 ± 0.17 

21.10 ± 0.29 
4.07 ± 0.13 
7.70 ± 0.19 
8.40 ± 0.19 

24.80 ± 0.31 
5.20 ± 0.19 

43.20 ± 0.56 
5.60 ± 0.13 

18.70 ± 0.04 
2.60 ± 0.13 
3.50 ± 0.13 

19.90 ± 0.38 

23.10 ± 0.69 

20.90 ± 1.13 
2.70 ± 0.13 
7.70 ± 0.13 

71.00 ± 1.25 

38.80 ± 0.50 

3.80 ± 0.13 
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: Sakamoto et al. (2008) 
Sakamoto et al. (2008) 



Golenetskii et al. (2005b) 

Golcnctskii et al. (200 [ 5a): 

Sakamoto et al. (2008) 

Krimm et al. (2009) : Sakamoto et al. (2008) 
Palmer et al. (2006) ; Sakamoto et al. (2008) 
Krimm et al. (2009) ; Sakamoto et al. (2008) 
Golenetskii et al. (2006d) ; Sakamoto et al. (2008) 
Krimm et al. (2009) ; Sakamoto et al. (2008) 
Krimm et al. (2009) ; Sakamoto et al. (2008) 
Sakamoto et al. (20081 

Golenetskii et al. (2006c1 : Sakamoto et al. (20081 
Golenetskii et al. (2006a1 ; Sakamoto et al. Q2008) 
Golenetskii et al. (2006b) ; Sakamoto et al. (2008) 
Sakamoto et al. (2008) 

Golenetskii et al. (2007b); Markwardt et al. (2007) 



Golenetskii et al. (2007a) ; Holland et al. (2007 ) 
Golenetskii et al. (2008f) ; Racusin et al. (2008) 
Golenetskii et al. (2008c) ; Stamatikos et al. (2008b) 
Golenetskii et al. (2008e) ; Sato et al. (2008) 
Krimm et al. (2009) ; Marshall et al. (2008) 
Krimm et al. (2009) ; Barthelmv et al. (20081 
Guidorzi et al. (20081 
Mangano et al. (20081 

Golenetskii et al. (2008d) ; Sbarufatti et al. (2008) 
Golenetskii et al. (2008a) ; Mangano et al. (2008) 
Golenetskii et al. (2008b) ; Marshall et al. (2008) 
Krimm et al. (20091 ; Ziaeepour et al. (20081 
Bissaldi &; McBreen (20081 ; Grupe et al. (2009) 
Connaughton (2009) ; Cannizzo et al. (2009 ) 
McBreen et al. (2009) ; Schadv et al. (2009 ) 
Golenetskii et al. (2009 ); Vetere et al. (2009) 



References — (1) IWatson et al.l 120061) : (2 ) IBerger fe Becked 120051) ; ( 3) IPiranomonte et all 120081): (4) IPenprase et al. | l|2006l); (5) 
IFvnbo et all 1120091); (6) IFvnbo et al.l 120091); (7) IProchaska et all ^20091: (8) IFvnbo et all 120091) :(9) IFvnbo et al.l 12009); (10) iFvnbo et"ahl 
1120091; (TTTlFvnbo e t aU 120091) ; (12) IFvnbo et al.l 1200911 : (13) IFvnbo et al.l 120091) ; ( 14) IJaunsen et al.l 120081) ; (15) ICenko et al.l 120071) : 
(16) Jakobsson ct al. (2007); (17) D'Elia et al. (2009); (18) Fynbo ct al. (2009); (19) Fynbo et al. (2009): (20) Fynbo et al. (2009); (21) 
Fynbo et al. (2009); (22) Cucchiara & Fox (2008); (23) Fynbo et al. (2009); (24) Fynbo ct al. (2009); (25) Prochaska ct al. (2009); (26) 
Fynbo et all I2OO9I); ( 27) IFvnbo et al.l 12009 ); (28) ICucchiara et al.l 120081) ; (29) Ichornock et al.l 12009) ; (30) ICenko et alj 1I2OO9I) ; (31) 
Wiersema et al.l 12009) . 

a l-second peak photon flux measured in photon cm -2 s -1 in the energy range 15 — 150 kcV. 

b Peak energy is given in keV. Values in brackets indicates estimated values using the method described in Sakamoto ct al. (2009). 

c Values in brackets indic ates estimated high-energy photon index, a, which is the mean value of the BATSE a distribution IKaneko et al.l 
120061 : [Sakamoto et aHl2009) . 

d Values in brackets indic ates estimated high-energy photon index, /3, which is the mean value of the BATSE (3 distribution IKaneko et al.l 
120061 : [Sakamoto et al]|2009) . 



Note. — Note that uncertainties of parameters that are reported with 90% confidence level have been reduced to lcr level for consistency. 
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Table 2 

grb redshift and calculated isotropic peak luminosity 



GRB 


Rodshift 


Peak Isotropic Luminosity a 


GRB 


Rodshift 


GRB050401 


2.899 


(1.3Btg;i!) 


X 


10 53 


GRB080319B 


0.937 


GRB050603 


2.821 


(6.32l° :«) 


X 


1Q 53 


GRB080319C 


1.949 


GRB050922C 


2.199 


l°" 1 '-0.01 1 


X 


1Q 52 


GRB080411 


1.030 


GRB051111 


1.550 




X 


1Q 52 


GRB080413A 


2.433 


GRB060206 


4.056 


(6.28+||) 


X 


10 52 


GRB080413B 


1.101 


GRB060210 


3.913 


(8.53_ g2 , 


X 


10 62 


GRB080430 


0.767 


GRB060418 


1.490 


(1.96±° ;f 3 ) 


X 


10 52 


GRB080603B 


2.689 


GRB060904B 


0.703 


(2.i8lS:SS) 


X 


10 61 


GRB080605 


1.640 


GRB060908 


1.884 


(1.54™) 


X 


10 52 


GRB080607 


3.036 


GRB060927 


5.464 


y-- 1 ' -o.io> 


X 


10 53 


GRB080721 


2.591 


GRB061007 


1.262 


(i-oii° : 2 o° 8 ) 


X 


10 53 


GRB080916A 


0.689 


GRB061021 


0.346 


(i.30±S;») 


X 


10 51 


GRB081222 


2.770 


GRB061121 


1.315 


(7.89± ;° 4 2 7 ) 


X 


10 62 


GRB090424 


0.544 


GRB070306 


1.496 


(8.67± 3 27 °) 


X 


10 51 


GRB090618 


0.540 


GRB071010B 


0.947 


(4.341;;") 


X 


10 51 


GRB090715B 


3.000 


GRB071020 


2.145 


(i.27iS:!S) 


X 


10 53 







Peak Isotropic Luminosity ' 



6.96^ 



+ 0.32s 
0.14/ 

6 04+ 8 04 ) 
5 49+ 1 11 ) 

5 38 +4 - 69 ~l 

- °°-0.S3> 

1 51+ 016l | 

1 - ° 1 -o.os> 

1 03+ 1 ' 30 "! 

2 QQ + 125s l 
, hk+0.56% 
1 - lo -0.09> 

' ■- Ly _0.41«' 

5 18+ 0S3 ~1 
o n + 19.90-, 

■ ,3U -0.15 > 

1 26 +0 07 l 
1 62+ 05 ) 
8 47+ 1 17 ) 

°- f " -0.34/ 

6 79+ 2 42 ) 



x 10 52 
x 10 52 
x 10 52 
x 10 52 
x 10 52 
x 10 51 
x 10 52 
x 10 53 
x 10 53 
x 10 53 
x 10 51 
x 10 53 
x 10 52 
x 10 51 
x 10 52 



a Isotropic equivalent peak photon luminosity in erg s 1 between GRB rest frame energy range 1 and 10,000 keV as described 
in Section [2] 




Fig. 1— The effect of noise on the CCF. Pan- 
els on the left show two synthetic light curves, in 
which a 10-second artificial lag is added. From 
top to bottom the noise level is increased as 0%, 
20%, 40% and 60% respectively. The correspond- 
ing CCF vs time delay plots are shown in the right 
panels along with Gaussian fits. 



CCFMax vs Noise Leve 



Noise Level (%) 
Lag vs Noise Leve 




Fig. 2. — The effect of noise on the maximum cor- 
relation of the CCF (CCFMax) and the extracted 
spectral lag. The horizontal line (blue) in the bot- 
tom panel indicates the 10-second artificial lag. 
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Table 3 

Spectral lag values of long duration Swift BAT GRBs 



GRB 


Iriggcr 1L) 


LagXX 


1 -f- A-S {S) 


T _1_ Y,. (v\ 
1 -f- AE {S) 


oin oizc ^msj 


T ^ 


( s ) 


LE (s) 


Lag Value fmsj 


GRB050401 


113120 


Lag 


21 


23.03 


29.43 


32 


-1 


.00 


1.50 


275 ± 131 






Lag 


31 


23.03 


29.43 


32 


-1 


.00 


2.00 


504 ± 117 






Lag 


41 


23.03 


29.43 


64 


-1 


.00 


2.00 


562 ± 140 






Lag 


32 


23.03 


29.43 


16 


-1 


.00 


1.00 


136 ± 87 






Lag 


42 


23.03 


29.43 


64 


-1 


.50 


1.50 


250 ± 112 






Lag 


43 


23.03 


29.43 


64 


-2 


.00 


2.00 


106 ±118 


GRB050603 


131560 


Lag 


21 


-3.83 


3.08 


8 


-0 


.40 


0.40 


46 ± 24 






Lag 


31 


-3.83 


3.08 


8 


-0 


.40 


0.40 


59 ± 22 






Lag 


41 


-3.83 


3.08 


16 


-0 


.40 


0.40 


86 ± 29 






Lag 


32 


-3.83 


3.08 


4 


-0 


.20 


0.20 


4 ± 11 






Lag 


42 


-3.83 


3.08 


16 


-0 


.40 


0.40 


34 ± 19 






Lag 


43 


-3.83 


3.08 


16 


-0 


.50 


0.50 


20 ± 18 


GRB050922C 


156467 


Lag 


21 


-2.70 


2.94 


8 


-0 


.40 


0.40 


9 ± 35 






Lag 


31 


-2.70 


2.94 


8 


-1 


.00 


1.00 


180 ± 50 






Lag 


41 


-2.70 


2.94 


16 


-1 


.00 


1.00 


188 ± 78 






Lag 


32 


-2.70 


2.94 


4 


-1 


.00 


1.00 


188 ± 39 






Lag 


42 


-2.70 


2.94 


16 


-1 


.00 


1.00 


178 ± 70 






Lag 


43 


-2.70 


2.94 


16 


-1 


.00 


1.00 


19 ± 72 


GRB051111 


163438 


Lag 


21 


-6.96 


28.62 


32 


-5 


.00 


4.00 


583 ± 273 






Lag 


31 


-6.96 


28.62 


32 


-4 


.00 


4.00 


1383 ± 288 






Lag 


41 


-6.96 


28.62 


128 


-4 


00 


8.00 


2343 ± 397 






Lag 


32 


-6.96 


28.62 


16 


-5 


.00 


4.00 


776 ± 200 






Lag 


42 


-6.96 


28.62 


64 


-5 


.00 


5.00 


1486 ± 314 






Lag 


43 


-6.96 


28.62 


64 


-5 


.00 


5.00 


866 ± 319 


GRB060206 


180455 


Lag 


21 


-1.29 


8.18 


8 


-1 


.50 


1.50 


241 ± 78 






Lag 


31 


-1.29 


8.18 


16 


-1 


.00 


2.00 


517 ± 85 






Lag 


41 


-1.29 


8.18 


64 


-1 


.50 


2.00 


331 ± 219 






Lag 


32 


-1.29 


8.18 


8 


-1 


.50 


1.50 


278 ± 74 






Lag 


42 


-1.29 


8.18 


64 


-1 


.50 


1.50 


82 ± 193 






Lag 


43 


-1.29 


8.18 


64 


-2 


.00 


2.00 


-163 ± 189 


GRB060210 


180977 


Lag 


21 


-3.37 


5.08 


64 


-5 


.00 


4.00 


700 ± 270 






Lag 


31 


-3.37 


5.08 


64 


-5 


.00 


4.00 


508 ± 254 






Lag 


41 


-3.37 


5.08 


256 


-5 


.00 


4.00 


1038 ± 324 






Lag 


32 


-3.37 


5.08 


64 


-4 


.00 


4.00 


-175 ± 174 






Lag 


42 


-3.37 


5.08 


128 


-4 


.00 


4.00 


98 ± 225 






Lag 


43 


-3.37 


5.08 


256 


-5 


.00 


2.00 


34 ± 195 


GRB060418 


205851 


Lag 


21 


-7.66 


33.04 


16 


-2 


.00 


2.00 


22 ± 62 






Lag 


31 


-7.66 


33.04 


32 


-2 


.00 


2.00 


109 ± 62 






Lag 


41 


-7.66 


33.04 


128 


-2 


.00 


2.00 


476 ± 196 






Lag 


32 


-7.66 


33.04 


16 


-2 


.00 


2.00 


87 ± 50 






Lag 


42 


-7.66 


33.04 


64 


-1 


.00 


1.00 


212 ± 100 






Lag 


43 


-7.66 


33.04 


64 


-1 


.00 


1.00 


162 ± 101 


GRB060904B 


228006 


Lag 


21 


-1.97 


10.32 


32 


-2 


.00 


2.00 


412 ± 195 






Lag 


31 


-1.97 


10.32 


32 


-2 


.00 


2.00 


560 ± 164 






Lag 


41 


-1.97 


10.32 


128 


-3 


.00 


3.00 


602 ± 296 






Lag 


32 


-1.97 


10.32 


32 


-2 


.00 


2.00 


247 ± 140 






Lag 


42 


-1.97 


10.32 


128 


-2 


.50 


3.00 


175 ± 292 






Lag 


43 


-1.97 


10.32 


128 


-3 


.00 


3.00 


32 ± 273 


GRB060908 


228581 


Lag 


21 


-10.91 


3.68 


32 


-2 


.00 


2.00 


118 ± 142 






Lag 


31 


-10.91 


3.68 


32 


-2 


.00 


2.00 


346 ± 185 






Lag 


41 


-10.91 


3.68 


128 


-4 


.00 


4.00 


367 ± 315 






Lag 


32 


-10.91 


3.68 


16 


-2 


.00 


2.00 


124 ± 86 






Lag 


42 


-10.91 


3.68 


64 


-2 


.00 


2.00 


233 ± 216 






Lag 


43 


-10.91 


3.68 


128 


-4 


.00 


4.00 


134 ± 253 


GRB060927 


231362 


Lag 


21 


-1.69 


8.04 


16 


-0 


.60 


0.60 


9 ±46 






Lag 


31 


-1.69 


8.04 


64 


-1 


.00 


1.00 


74 ± 62 






Lag 


41 


-1.69 


8.04 


256 


-1 


.50 


1.50 


200 ± 133 






Lag 


32 


-1.69 


8.04 


16 


-1 


.00 


1.00 


103 ± 45 






Lag 


42 


-1.69 


8.04 


128 


-1 


.20 


1.50 


229 ± 112 
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Fig. 3. — Swift-BAT prompt gamma-ray (8 ms 
time bin) light curves for GRB 060206 with canon- 
ical energy channels 2 (25-50 keV) and 3 (50-100 
keV). 



100 200 300 400 500 600 

Spectral Lag (ms) 

Fig. 5. — Histogram of 1,000 simulated spectral 
lag values. We take the standard deviation of the 
distribution of simulated spectral lag values as the 
uncertainty of the fitted spectral lag value which 
was found in Fig. [4] The final spectral lag value 
is 278 ± 74 ms. 



Gaussian Fit Lag : 278. ± 13. m 




Fig. 4. — CCF as a function of time delay for 
the two light curves in Fig|3j The time delay 
corresponding to the peak of the Gaussian fit is 
the spectral lag of the burst, which is 278 ± 13 
ms. The uncertainty quoted here is from the fit, 
which tends to be factor of two or more less than 
the value obtained by the Monte Carlo simulation 
shown in Fig. [5j 




1.0 1.5 2.0 2.5 

log Lag 21/(1 +z) msec 



Fig. 6. — Isotropic luminosity as a function of 
spectral lag between BAT channel 2 (25-50 keV) 
and 1 (15-25 keV). The "Gold", "Silver", and 
"Bronze" samples are represented with red circles, 
blue diamonds, and green triangles respectively. 
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-) 1 1 1 1 1 1 1 1 1 1 1 1 1- — 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



log Lag 32/(1 +z) msec 

Fig. 7. — Isotropic luminosity as a function of 
spectral lag between BAT channel 3 (50-100 keV) 
and 2 (25-50 keV). The "Gold", "Silver", and 
"Bronze" samples are represented with red circles, 
blue diamonds, and green triangles respectively. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

log Lag 43/(1 +z) msec 



Fig. 9. — Isotropic luminosity as a function of 
spectral lag between BAT channel 4 (100-200 keV) 
and 3 (50-100 keV). The "Gold", "Silver", and 
"Bronze" samples are represented with red circles, 
blue diamonds, and green triangles respectively. 
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50-| , 1 , 1 , 1 , 1 , 1 , 1 H 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

log Lag 31/(1+z) msec 

Fig. 8. — Isotropic luminosity as a function of 
spectral lag between BAT channel 3 (50-100 keV) 
and 1 (15-25 keV). The "Gold", "Silver", and 
"Bronze" samples are represented with red circles, 
blue diamonds, and green triangles respectively. 



'"• t- m i ■ 1 


1 1 1 1 1 1 




Mlll 'Toei222 

fife* 


090424 - 














%. 






0610^J \^ ^ 





0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

log Lag 42/(1 +z) msec 



Fig. 10. — Isotropic luminosity as a function of 
spectral lag between BAT channel 4 (100-200 keV) 
and 2 (25-50 keV). The "Gold", "Silver", and 
"Bronze" samples are represented with red circles, 
blue diamonds, and green triangles respectively. 
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-) 1 1 1 1 1 1 1 1 1 1 1 1 1- — 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



log Lag 41/(1+z) msec 

Fig. 11. — Isotropic luminosity as a function of 
spectral lag between BAT channel 4 (100-200 keV) 
and 1 (15-25 keV). The "Gold", "Silver", and 
"Bronze" samples are represented with red circles, 
blue diamonds, and green triangles respectively. 
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Table 3 — Continued 



GRB 


Trigger ID 


LagXX 


T + X S (s) 


T + X E (s) 


Bin Size (ms) 


LS (s) 


LE (s) 


Lag Value (ms' 






Lag 43 


-1.69 


8.04 


128 


-1 .20 


1.50 


i oc _l_ i n i 
Izu ± 1U1 




ZoZooo 


Lag 21 


23.86 


65.08 


o 
Z 


-0.30 


0.50 


101 i 17 






Lag 31 


23.86 


65.08 


2 


-0.30 


0.50 


154 i 19 






Lag 41 


23.86 


65.08 


4 


-0.50 


0.80 


HOC _|_ OO 

ZOO It ZO 






Lag 32 


23.86 


65.08 


2 


-0.20 


0.20 


30 i 8 






Lag 42 


23.86 


65.08 


2 


-0.40 


0.40 


129 i 17 






Lag 43 


23.86 


65.08 


2 


-0.30 


0.40 


82 i 9 


PDRflfil noi 




Lag 21 


-0.46 


14.64 


8 


-1.00 


1.00 


— 25 i 52 






Lag 31 


-0.46 


14.64 


o 


-1.00 


1.00 


49 i 51 






Lag 41 


-0.46 


14.64 


32 


-1.60 


1.60 


o on _l_ o c 
zoy ± 8o 






Lag 32 


-0.46 


14.64 


8 


-1.00 


1.00 


62 i 42 






Lag 42 


-0.46 


1 A P.A 


OZ 


1 HA 
-l.UU 


1.20 


OAH -4- la 
Z4o It to 






Lag 43 


n A k 

-U.40 


1 A R.A 

14.04 


qo 
OZ 


1 HA 
-l.UU 


1 on 
l.ZU 


loo it i y 


\j t\.r> U u llzl 


zoyoyy 


T it oi 

i^ag zi 


0U.44 


an p.p. 

OU.OO 


1 
1 


n on 
-U.zU 


n oa 
u.zu 


io J. iq 
lo it lo 






Lag 31 


60.44 


80.66 


1 


-0.20 


0.20 


16 it 12 






Lag 41 


60.44 


80.66 


4 


-0.40 


0.40 


26 i 26 






Lag 32 


60.44 


80.66 


1 


-0.20 


0.25 


17 i 7 






Lag 42 


60.44 


80.66 


2 


-0.20 


0.25 


28 i 12 






Lag 43 


60.44 


80.66 


z 


-0.20 


0.25 


25 ill 


T> Dn7flQflC 

^jHriU /UoUD 


zooobl 


Lag 21 


90.00 


1 18.42 


8 


-2.00 


2.00 


08 ± lUb 






Lag 31 


on nn 
yu.uu 


1 1 Q A O 


1 K 
ID 


O OA 
-z.UU 


o nn 
z.UU 


1 Ad -l- i nn 
14D It 1UU 






Lag 41 


on nn 
yu.uu 


1 1 Q A O 
1 10.4Z 


RA 
D4 


A nn 
-4.UU 


^ nn 
D.UU 


1 noc _l_ qoi 

luoo it oyi 






T it QO 

.Lag oz 


on nn 
yu.uu 


1 1 Q A O 
1 10.4Z 


Q 

o 


o nn 
-z.UU 


o nn 
z.UU 


i i a 4- i no 
1 14 ± lUz 






Lag 42 


90.00 


1 18.42 


64 


-4.00 


6.00 


i noc _l_ Qon 

iuy8 ± oyy 






Lag 43 


90.00 


1 18.42 


64 


-4.00 


6.00 


nnn _l_ Ana 
9UU ± 4Uo 


r'PRfi7i ni nR 
OrvxSU < 1U1U.D 


zyo i yo 


Lag 21 


-1.70 


17.24 


2 


-1.00 


1.00 


— 26 i 48 






Lag 31 


-1.70 


17.24 


4 


-1.00 


1.00 


146 i 52 






Lag 41 


-1.70 


17.24 


32 


-2.00 


4.00 


ino/i -U i tzi 

1UZ4 It IDO 






Lag 32 


-1.70 


17.24 


4 


-1.00 


1.00 


185 i 47 






Lag 42 


-1.70 


17.24 


32 


-2.00 


4.00 


i nn i" i it? 
1UU5 ± lot 






Lag 43 


-1.70 


17.24 


32 


-2.00 


4.00 


745 i 161 


r'RRfl71 non 
vjr\£3U / lUzU 




Lag 21 


-O.ZZ 


1. 14 


z 


n in 
-U. 1U 


n i ^ 

U. 10 


I It I 






T „ _. q 1 

ijag oi 


q no 
-O.ZZ 


11/1 
±. ±4 


z 


n in 
-U. 1U 


n on 
U.zU 


07 _l_ 1 o 
o / it Iz 






Lag 41 


-3.22 


1. 14 


8 


-0.50 


0.50 


— 50 i 30 






Lag 32 


-3.22 


1. 14 


2 


-0.10 


0.25 


47 i 7 






Lag 42 


q no 
-O.ZZ 


11/1 
1. 14 


4 


-0.10 


0.30 


69 i 12 






Lag 43 


-3.22 


11/1 
1. 14 


4 


-0.20 


0.30 


28 i 9 


GRB080319B 


306757 


Lag 21 


-2.85 


57.57 


2 


-0.10 


0.14 


15 it 2 






Lag 31 


-2.85 


57.57 


2 


-0.10 


0.14 


32 ±3 






Lag 41 


-2.85 


57.57 


2 


-0.20 


0.20 


80 ib 17 






Lag 32 


-2.85 


57.57 


2 


-0.10 


0.14 


23 ±2 






Lag 42 


-2.85 


57.57 


2 


-0.20 


0.30 


88 ±8 






Lag 43 


-2.85 


57.57 


2 


-0.20 


0.20 


26 ±5 


GRB080319C 


306778 


Lag 21 


-0.77 


13.31 


16 


-1.00 


1.00 


106 ± 78 






Lag 31 


-0.77 


13.31 


16 


-2.00 


2.00 


216 ± 70 






Lag 41 


-0.77 


13.31 


64 


-2.00 


2.00 


89 ib 132 






Lag 32 


-0.77 


13.31 


16 


-1.00 


1.00 


134 ± 58 






Lag 42 


-0.77 


13.31 


32 


-1.00 


1.00 


-77 ib 95 






Lag 43 


-0.77 


13.31 


32 


-1.00 


1.00 


-119 ±99 


GRB080411 


309010 


Lag 21 


38.46 


48.45 


2 


-1.00 


1.00 


103 ± 12 






Lag 31 


38.46 


48.45 


2 


-1.00 


1.00 


220 ± 13 






Lag 41 


38.46 


48.45 


2 


-1.00 


1.00 


322 ± 27 






Lag 32 


38.46 


48.45 


2 


-1.00 


1.00 


122 ± 11 






Lag 42 


38.46 


48.45 


2 


-1.00 


1.00 


230 ± 26 






Lag 43 


38.46 


48.45 


2 


-1.00 


1.00 


112 ± 26 


GRB080413A 


309096 


Lag 21 


-0.42 


9.05 


8 


-1.00 


1.00 


96 ± 60 






Lag 31 


-0.42 


9.05 


8 


-1.00 


1.00 


242 ±65 






Lag 41 


-0.42 


9.05 


64 


-1.00 


2.00 


542 ± 125 






Lag 32 


-0.42 


9.05 


8 


-1.00 


1.00 


157 ± 43 
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Table 3 — Continued 



GRB 


Trigger ID 


LagXX 


T + X S (s) 


T + X E (s) 


Bin Size (ms) 


LS (s) 


LE (s) 


Lag Value (ms 






Lag 42 


-0.42 


9.05 


32 


-1.00 


2.00 


418 i 111 






Lag 43 


-0.42 


9.05 


32 


-1 .00 


2.00 


oa o -I- ins 
z4y it ±Uo 


\j tx ±5 U o U 4 ± o ±3 


ouy i i i 


Lag 21 


1 A A 
-±.44 


4.96 


8 


-1.00 


1.00 


59 ± 35 






Lag 31 


1 A A 
-±.44 


4.96 


8 


-1.00 


1 .00 


1/1/14- q? 
±44 ± or 






Lag 41 


-1.44 


4.96 


16 


-1.00 


1.00 


q k q 4- fit; 
ooo ± 06 






Lag 32 


-1.44 


4.96 


8 


-1.00 


1.00 


82 ± 28 






T it A 

i^ag 4z 


1 A A 
-±.44 


4.96 


16 


-1.00 


1.00 


276 ± 56 






Lag 43 


1 A A 
-±.44 


4.96 


16 


-1.00 


1.00 


188 ±55 


ptj Rrififi 1 *?n 
vj rtrS U o U4o U 


q i r\f? i q 


Lag 21 


-1.24 


12.84 


32 


-2.00 


2.00 


T7n 4- Qft 
z / U It OO 






Lag 31 


-1.24 


12.84 


32 


-2.00 


2.00 


oni _l_ i no 
oyl ± luy 






Lag 41 


-1.24 


12.84 


256 


-4.00 


4.00 


Tqn _l_ ota 
( oil ± o i 4 






Lag 32 


-1.24 


12.84 


32 


-2.00 


2.00 


oo i 1 nn 
oo ± 1UU 






Lag 42 


1 OA 


±Z.04 


ZOK) 


A nn 
-4.UU 


A nn 
4.UU 


04U It Oo I 






ijag 4o 


1 OA 
-±.Z4 


1 O S/f 
±Z.04 


O^K 

ZOO 


A nn 
-4.UU 


A nn 
4.UU 


qoo 4_ on? 
ooo It oy I 




oloOo ( 


Lag 21 


-0.54 


5.10 


16 


-1.00 


1.00 


noi 4_ p. i 
— zzz it 0± 






Lag 31 


-0.54 


5.10 


16 


-1.00 


1.00 


1 07 4- C7 

— ly ( ± o ( 






Lag 41 


-0.54 


5.10 


32 


-1.00 


1.00 


A oi 4_ i cq 
— 4Z i it ±0o 






Lag 32 


-0.54 


5.10 


16 


-1.00 


1.00 


50 ± 41 






Lag 42 


n ka 

-U.04 


5.10 


oZ 


-1.00 


0.50 


— 103 ± 71 






Lag 43 


-0.54 


5.10 


32 


-1.00 


0.50 


— 172 ± 56 


*-jrtr>UoUDUO 




Lag 21 


-0.40 


lo. oo 


A 

4 


i nn 


i nn 


4- oo 
Oo it zy 






Lag 31 


-0.40 


10. Oo 


A 

4 


i nn 
-±.uu 


i nn 
±.UU 


qq _|_ qq 

yo it oo 






T „„ 11 

.Lag 4± 


-0.40 


10. Oo 


±0 


n ^n 
-U.OU 


i on 


1 4_ on 
±yo it oy 






Lag 32 


-5.46 


15.53 


2 


-0.30 


0.40 


73 ± 11 






Lag 42 


-5.46 


15.53 


8 


-0.30 


0.40 


96 ± 17 






Lag 43 


-0.40 


15.53 


8 


-0.30 


0.40 


39 ± 12 


UnljUoUDU / 


313417 


Lag 21 


-6.13 


12.05 


8 


-0.40 


0.40 


121 ± 119 






Lag 31 


-6.13 


1 o nr^ 


8 


-0.40 


0.60 


1 4- ICk 

±oo it oy 






Lag 41 


-6.13 


12.05 


16 


-0.40 


0.60 


194 ± 43 






Lag 32 


-6.13 


12.05 


8 


-0.40 


0.40 


19 ± 17 






Lag 42 


-6.13 


12.05 


8 


-0.40 


0.40 


64 ± 23 






Lag 43 


-o. ±o 




Q 
O 


n a n 
-U.4U 


n a n 
U.4U 


Or^ 4_ i q 
ZO It ±o 




ol/ OUo 


T it oi 

ijag zi 


q qn 

-o.oy 


a p.a 
0.04 


RA 


o nn 
-Z.Vl) 


o nn 
z.UU 


QQ 4- 1 AQ 

yy it ±4y 






Lag 31 


-3.39 


8.64 


64 


-2.00 


2.00 


1 OO 1 1 no 
LZZ ± loo 






Lag 41 


-3.39 


8.64 


128 


-2.00 


2.00 


q a i 4_ i on 
o4± it loz 






Lag 32 


-3.39 


8.64 


16 


-0.80 


0.80 


16 ± 58 






Lag 42 


-3.39 


8.64 


32 


-0.80 


0.80 


O ^ft 4- QflQ 
ZOO it OUO 






Lag 43 


-3.39 


8.64 


32 


-0.80 


0.80 


167 ± 69 


GRB080916A 


324895 


Lag 21 


-2.66 


39.58 


16 


-2.00 


3.00 


566 ± 172 






Lag 31 


-2.66 


39.58 


32 


-2.00 


4.00 


1468 ± 202 






Lag 41 


-2.66 


39.58 


256 


-4.00 


6.00 


2879 ± 271 






Lag 32 


-2.66 


39.58 


32 


-4.00 


4.00 


821 ± 100 






Lag 42 


-2.66 


39.58 


128 


-2.00 


6.00 


1900 ± 165 






Lag 43 


-2.66 


39.58 


64 


-2.00 


4.00 


842 ± 143 


GRB081222 


337914 


Lag 21 


-0.80 


15.58 


2 


-0.80 


1.20 


127 ± 41 






Lag 31 


-0.80 


15.58 


2 


-0.80 


1.00 


262 ± 47 






Lag 41 


-0.80 


15.58 


16 


-2.00 


3.00 


610 ± 111 






Lag 32 


-0.80 


15.58 


2 


-0.80 


0.80 


113 ±30 






Lag 42 


-0.80 


15.58 


8 


-1.80 


2.80 


444 ± 107 






Lag 43 


-0.80 


15.58 


8 


-1.80 


2.80 


197 ± 110 


GRB090424 


350311 


Lag 21 


-0.94 


4.95 


1 


-0.10 


0.25 


20 ± 12 






Lag 31 


-0.94 


4.95 


2 


-0.10 


0.25 


29 ± 13 






Lag 41 


-0.94 


4.95 


4 


-0.10 


0.25 


39 ± 15 






Lag 32 


-0.94 


4.95 


1 


-0.10 


0.25 


23 ±9 






Lag 42 


-0.94 


4.95 


4 


-0.10 


0.25 


27 ± 13 






Lag 43 


-0.94 


4.95 


4 


-0.20 


0.30 


17 ±9 


GRB090618 


355083 


Lag 21 


46.01 


135.35 


4 


-1.00 


1.00 


255 ±21 






Lag 31 


46.01 


135.35 


4 


-1.00 


1.00 


447 ± 26 






Lag 41 


46.01 


135.35 


4 


-1.00 


2.00 


894 ± 43 
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Table 3 — Continued 



GRB Trigger ID 


LagXX 


T + X S (s) 


T + X E (s) 


Bin Size (ms) 


LS (s) 


LE (s) 


Lag Value (ms) 




Lag 32 


46.01 


135.35 


4 


-1.00 


1.00 


173 ± 18 




Lag 42 


46.01 


135.35 


4 


-1.00 


1.50 


483 ± 34 




Lag 43 


46.01 


135.35 


4 


-1.00 


1.50 


283 ± 34 


GRB090715B 357512 


Lag 21 


-4.80 


21.06 


16 


-2.50 


2.50 


288 ± 117 




Lag 31 


-4.80 


21.06 


16 


-2.50 


2.50 


732 ± 127 




Lag 41 


-4.80 


21.06 


64 


-2.50 


3.00 


1080 ± 224 




Lag 32 


-4.80 


21.06 


8 


-2.50 


2.50 


470 ± 100 




Lag 42 


-4.80 


21.06 


32 


-2.50 


2.50 


928 ± 229 




Lag 43 


-4.80 


21.06 


32 


-2.50 


2.50 


375 ± 215 



Table 4 

Correlation coefficients and fit parameters 



Channels 


Correlation 






Best Fit 






X 2 /ndf 


Channel 21 


-0.63 


log L iso 


= (54.8 ± 0.2) 


- (1.4 ± 0.1) log 


La S 21(l + zY 


l 


189.4/19 


Channel 32 


-0.66 


log L iso 


= (54.5 ± 0.2) 


- (1.2 ± 0.1) log 


Lag32(l + z)~ 


l 


216/25 


Channel 31 


-0.60 


log L iso 


= (55.5 ± 0.2) 


- (1.5 ± 0.1) log 


La 9 31(l + z)- 


l 


410.8/26 


Channel 43 


-0.77 


log L iso 


= (55.0 ± 0.3) 


- (1.4 ± 0.1) log 


La S 43(l + z)- 


l 


109/20 


Channel 42 


-0.75 


log L iso 


= (55.4 ± 0.1) 


- (1.4 ± 0.1) logLa 9 42(l + z)~ 


l 


178.8/23 


Channel 41 


-0.67 


log L iao 


= (56.7 ± 0.3) 


- (1.8 ± 0.1) log 


La 9 41(l + z) 


l 


212.1/25 



Table 5 

Correlation coefficients with various corrections 



Channels 


No correction 


z-correction 


k-correction 


both corrections 


Channel 21 


-0.38 


-0.63 


-0.29 


-0.55 


Channel 32 


-0.43 


-0.66 


-0.33 


-0.60 


Channel 31 


-0.39 


-0.60 


-0.31 


-0.54 


Channel 43 


-0.61 


-0.77 


-0.54 


-0.73 


Channel 42 


-0.58 


-0.75 


-0.51 


-0.71 


Channel 41 


-0.43 


-0.67 


-0.32 


-0.61 
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